home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / libblas / spr.z / spr
Encoding:
Text File  |  2002-10-03  |  3.4 KB  |  85 lines
  1. SPR(3F)                                               Last changed: 11-2-98
  2.  
  3.  
  4. NNAAMMEE
  5.      SSSSPPRR, DDSSPPRR - Performs symmetric rank 1 update of a real or complex
  6.      symmetric packed matrix
  7.  
  8. SSYYNNOOPPSSIISS
  9.      Real
  10.  
  11.         CCAALLLL SSSSPPRR ((_u_p_l_o,, _n,, _a_l_p_h_a,, _x,, _i_n_c_x,, _a_p))
  12.  
  13.      Double precision
  14.  
  15.         CCAALLLL DDSSPPRR ((_u_p_l_o,, _n,, _a_l_p_h_a,, _x,, _i_n_c_x,, _a_p))
  16.  
  17. IIMMPPLLEEMMEENNTTAATTIIOONN
  18.      IRIX systems
  19.  
  20. DDEESSCCRRIIPPTTIIOONN
  21.      These routines perform the following symmetric rank 1 operation:
  22.                        _T
  23.           _A <- _a_l_p_h_a _x_x  + _A
  24.             _T
  25.      where _x   is the transpose of _x, _a_l_p_h_a is a real or complex scalar, _x
  26.      is an _n-element vector, and _A is an _n-by-_n symmetric packed matrix.
  27.  
  28.      These routines have the following arguments:
  29.  
  30.      _u_p_l_o      Character*1.  (input)
  31.                Specifies whether the upper or lower triangular part of
  32.                matrix _A is packed into the array argument _a_p, as follows:
  33.  
  34.                _u_p_l_o= 'U' or 'u':  the upper triangular part of _A is being
  35.                supplied in the argument _a_p.
  36.                _u_p_l_o= 'L' or 'l':  the lower triangular part of _A is being
  37.                supplied in the argument _a_p.
  38.  
  39.      _n         Integer.  (input)
  40.                Specifies the order of matrix _A.  _n >= 0.
  41.  
  42.      _a_l_p_h_a     Scalar alpha.  (input)
  43.                SSSSPPRR: Real.
  44.                DDSSPPRR: Double precision.
  45.  
  46.      _x         Array of dimension 1+(_n-1) * |_i_n_c_x|.  (input)
  47.                SSSSPPRR: Real array.
  48.                DDSSPPRR: Double precision array.
  49.                Contains vector _x.
  50.  
  51.      _i_n_c_x      Integer.  (input)
  52.                Specifies the increment for the elements of _x.  _i_n_c_x must
  53.                not be 0.
  54.  
  55.      _a_p        Array of dimension (_n(_n+1))/2 .  (input and output)
  56.                SSSSPPRR: Real array.
  57.                DDSSPPRR: Double precision array.
  58.  
  59.                Before entry with _u_p_l_o = 'U' or 'u', array _a_p must contain
  60.                the upper triangular part of the symmetric matrix packed
  61.                sequentially, column-by-column, so that _a_p(1) contains
  62.                _A(1,1), _a_p(2) contains _A(1,2), _a_p(3) contains _A(2,2), and so
  63.                on.  On exit, the upper triangular part of the updated
  64.                matrix overwrites array _a_p.
  65.  
  66.                Before entry with _u_p_l_o = 'L' or 'l', array _a_p must contain
  67.                the lower triangular part of the symmetric matrix packed
  68.                sequentially, column-by-column, so that _a_p(1) contains
  69.                _A(1,1), _a_p(2) contains _A(2,1), _a_p(3) contains _A(3,1), and so
  70.                on.  On exit, the lower triangular part of the updated
  71.                matrix overwrites array _a_p.
  72.  
  73. NNOOTTEESS
  74.      SSSSPPRR/DDSSPPRR is a Level 2 Basic Linear Algebra Subprogram (Level 2 BLAS).
  75.  
  76.      When working backward (_i_n_c_x < 0), each routine starts at the end of
  77.      the vector and moves backward, as follows:
  78.  
  79.           _x(1-_i_n_c_x * (_n-1)), _x(1-_i_n_c_x * (_n-2)), ..., _x(1)
  80.  
  81. SSEEEE AALLSSOO
  82.      HHPPRR(3F)
  83.  
  84.      This man page is available only online.
  85.